Mapping the Evolution of Optically-Generated Rotational Wavepackets 
in a Room Temperature Ensemble of D 2 
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A coherent superposition of rotational states in D2 has been excited by nonresonant ultrafast 
(12 femtosecond) intense (2 x 10 14 Wcm -2 ) 800 nm laser pulses leading to impulsive dynamic 
alignment. Field-free evolution of this rotational wavepacket has been mapped to high temporal 
resolution by a time-delayed pulse, initiating rapid double ionization, which is highly sensitive to 
the angle of orientation of the molecular axis with respect to the polarization direction, 9. The 
detailed fractional revivals of the neutral D2 wavepacket as a function of 9 and evolution time have 
been observed and modelled theoretically. 

PACS numbers: 42.50.Hz, 33.80.Gj, 33.80.Wz 
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The importance of being able to enforce spatial order 
on an initially random ensemble of molecules has been 
recognized since the discovery of steric effects in chemi- 
cal reactions [l[ . Traditional 'brute force' techniques em- 
ploying strong DC fields 0, HJ have recently yielded to 
new, more subtle and yet more powerful and versatile 
techniques using intense laser systems 0, S| • In partic- 
ular intense femtosecond pulses have been successfully 
used to align molecules along an axis for application 
in areas such as the study of fragmentation dynamics 
i 0, high harmonic generation | @ and the creation 
of single cycle laser pulses fioj . In this case the align- 
ment of a molecule (or spatial ordering of an ensemble 
of molecules) evolves temporally following the creation 
of a rotational wavepacket by a linearly polarized laser 
pulse of far shorter duration than the natural period of 
rotation. 

In this Letter, we report on a detailed ultrafast study 
of the temporal evolution of such a rotational wavepacket 
in the neutral deuterium molecule D2. The wavepacket 
is created impulsively by a 12 femtosecond laser pulse, 
with temporal and angular ordering probed at some later 
time with a similar pulse that initiates sequential double 
ionization (SI) i.e. D 2 -> -> D+ + D+ [11]. Such 
high resolution measurements represent the state-of-the- 
art in ultrafast molecular physics in the high rotational 
frequency limit as dictated by this most fundamental and 
theoretically tractable of molecules. 

It is known that an intense linearly polarized laser 
pulse interacting with an ensemble of molecules gener- 
ates a degree of spatial alignment from a random ensem- 
ble 0, [B| . The only condition for this phenomenon is that 
the molecular polarizability is anisotropic. Thus, the in- 
duced dipole moment created in the interaction with the 
electric field of the laser generates a torque that causes 
the axis of maximum polarizability to librate around the 



polarization vector of the laser field. Quantum mechani- 
cally, the process is described as a sequence of Rabi-type 
cycles accompanied by the exchange of two quanta of 
angular momentum between the molecule and the laser 
field. If the aligning pulse duration is much greater than 
the rotational period, the system evolves adiabatically 
to the the so-called pendular states, which dissipate as 
soon as the laser pulse passes [12] • However, when the 
pulse duration is much shorter than the rotational pe- 
riod, a superposition of rotational states is created that 
outlives the laser pulse, and the phases in the rotational 
wavepacket continue to evolve [5J. As a consequence of 
the quantization of the J states, a periodic dephasing 
and rephasing of the wavepacket occurs in the field-free 
regime. At integer and half- integer multiples of (2Bc)~ 1 
(where B is the rotational constant and c is the speed of 
light in vacuum) a revival of the rotational wavepacket is 
expected, at which point the ensemble exhibits a signifi- 
cant degree of alignment, with the molecular axes parallel 
to the polarization vector of the initial pulse. The rela- 
tive phases of the J states also produces antialignment in 
the ensemble, whereby the molecular axes are orientated 
to lie on or near the plane normal to the initial pulse po- 
larization. Alignment and antialignment are well distin- 
guished temporally. As a consequence of the nuclear spin 
statistics, homonuclear systems can also present partial 
revivals at 1/4 and 3/4 of the rotational period. 

Previous observations of impulsive alignment have 
been restricted by the non-availability of laser pulses suf- 
ficiently shorter than the rotational periods. However, 
rotational wavepackets have been previously observed in 
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171 ] etc.), where rotations occur with 
periods on the picosecond time scale. The rotational re- 
vival in D 2 occurs at (2Bc)^ 1 = 558 fs, requiring an align- 
ing pulse of the order of tens of femtoseconds in duration. 
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FIG. 1: The experimental configuration used to observe rota- 
tional wavepackets in D2 imprinted on the D + rapid sequen- 
tial ionization (SI) yield (D2 — > — > D + + D + ) as measured 
with a time-of-flight mass spectrometer (TOFMS). The time 
delay between the pump and probe pulses is At, and the angle 
9 between pump and probe polarization directions is varied 
by rotating the pump polarization. As shown, 6 — 0. 



Dynamic alignment of D2 using 10 fs pulses has recently 
been demonstrated by Lee et al [18J , however, these ob- 
servations were of limited temporal resolution and the 
angular evolution of the rotational wavepacket was not 
measured. 

In the present work, the spatio-temporal evolution of 
the D2 rotational wavepacket has been simulated follow- 
ing the procedure described in (l9j . Briefly, a thermal en- 
semble of rigid rotors is considered, each of which gives 
rise to a superposition of rotational states, vP^m,) = 
Sj>|Mj| m which Fj i j{t) are the time- 

dependent complex coefficients (note that M; is con- 
served in a linearly polarized field). These coefficients 
are calculated by numerically solving the time-dependent 
Schrodinger equation, and are used to calculate the ther- 
mal average of the degree of alignment and the angular 
evolution of the ensemble. 

Due to the bosonic character of D2, the total wavefunc- 
tion of the molecule must be symmetric under inversion. 
The electronic ground state is symmetric, and the nuclei 
can form six symmetric (ortho) and three antisymmet- 
ric (para) nuclear wavefunctions. This restricts the J 
states populated: para-D 2 occupies only odd J-states, 
and ortho-D 2 occupies only even J-states. This results 
in a 2 : 1 weighting (ortho-D2 : para-D2) in a thermal 
ensemble. For a room temperature sample we expect the 
following initial populations (in parenthesis) for the dif- 
ferent J states: J = (0.185), 1 (0.208), 2 (0.386), 3 
(0.112), 4 (0.0899), 5 (0.0128) and 6 (0.00522). 

The experiment was carried out using 800nm, 0.4 mJ, 
10 fs pulses at 1 kHz repetition rate. The pulses were 
split by a 4 fan thick pellicle beamsplitter. A co-linear 




FIG. 2: (a) False colour representation of the D + rapid se- 
quential ionization (SI) yield (D2 — * — > D + + D + ) as 
a function of pump-probe delay (thus wavepacket evolution 
time) At and ion kinetic energy, (b) Integrated SI yield com- 
pared to theoretically predicted influence of an impulsively 
excited rotational wavepacket in D2 (red). Throughout, 8 — 
7r/2, see Fig. 1, and the step in At — 1 fs. (c) Integrated 
SI yield following an 11-point smoothing. The 1/4, 1/2 and 
3/4 revivals are more clearly observed. Fourier Transform 
of the theoretical simulation and experimental SI yield are 
shown in (d) and (e) respectively. The frequency components 
present in the rotational wavepacket are beats between rota- 
tional states AJ and AJ+2. The Fourier amplitudes thus 
return a measure of rotational population. 



Mach-Zehnder interferometer configuration was used to 
delay the resulting 70 /iJ, 12 fs probe pulse with respect 
to the 73_uJ, 12 fs pump pulse with 300 attosecond res- 
olution [20j. A half- wave plate mounted in one arm of 
the interferometer allowed the pump polarization to be 
rotated through an angle 9 with respect to the probe 
polarization, which was fixed parallel to the spectrome- 
ter axis. In comparison to the core-annulus method [ljj ]. 
the present technique does suffer greater energy loss in 
the split pulses. However, this is compensated by the 
temporal consistency of the spatial overlap at focus as a 
function of time delay over a far greater range. Further- 
more the present technique allows polarization control. 
These features are essential to the present study. 

Following transmission through a fused-silica window 
into an ultra-high vacuum chamber, the pump-probe 
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pulses are reflection focused with a silver-coated spherical 
mirror capable of supporting the bandwidth of our laser 
pulses without introducing group-velocity dispersion, as 
shown in Fig. 1. The laser focus is situated in the source 
region of a tightly apertured (250 fim) ion time-of-flight 
mass spectrometer (TOFMS) [20j]. The aperture serves 
two purposes: to limit the angular acceptance (< 2 de- 
grees for a 5.5 eV D + ion), and to dramatically limit the 
field of view of the spectrometer. By observing only the 
central 250 /an 'slice' of the f/5 focus, the instrument is 
sensitive only to those molecules exposed to a narrow in- 
tensity range. As a consequence, spatial integration over 
the focal volume is obviated, and the behaviour of the 
rotational wavepacket is clearly observed. This results in 
a remarkably high resolution measurement of the rota- 
tional wavepacket as compared to earlier observations. 

Fig. 2(a) shows the measured rapid sequential ioniza- 
tion D + yield in false colour as a function of pump-probe 
delay (At, also referred to as evolution time) and D + ion 
kinetic energy; here, 9 = ir/2. We are interested in the 
modulation in the SI yield < At < 800 fs, the signa- 
ture of the revival of a rotational wavepacket in the D2 
molecule. Importantly, in Fig. 2, the orthogonality of 
the aligning pump pulse to the TOFMS (and probe po- 
larization) axis results in the initial maximal alignment 
of the ensemble producing a minimum in the SI signal as 
the majority of the molecules are aligned perpendicular 
to the detector. 

The solid circles in Fig. 2(b) represent the exper- 
imentally obtained high kinetic energy D + yield from 
the rapid SI of D2 at small internuclear separation for 
9 = tt/2. The half- (280 fs) and full-revivals (560 fs) 
are clearly apparent in the results of both the integrated 
SI yield and the theoretical calculation, with excellent 
agreement observed between experiment and theory. At 
small delays < At < 100 fs, the discrepancy between 
experiment and theory is the result of the temporal wings 
of the pump and probe pulses overlapping to produce el- 
liptical/circular polarization depending on At. When At 
> 200fs, the time-dependent structure of the SI yield is 
purely the result of the rotational wavepacket. 

In Fig. 2(b), between 200 < At < 800 fs, there is 
evidence of fine structure in the experimentally observed 
wavepacket dependence, the result of the 1/4, 1/2 and 
3/4 revivals at respective fractions of (25c) _1 . Since the 
D + yield is recorded every 1 fs, an 11-point smoothing 
algorithm applied to the integrated SI yield (Fig. 2b) 
suppresses statistical fluctuations and hence enhances the 
fidelity of the fine structure observation without loosing 
temporal resolution, as shown in Fig. 2(c). 

Through the Fourier transform (FT) of the experimen- 
tal signal (Fig. 2b), we recover the frequency components 
present in the rotational wavepacket. Fig. 2(d) and Fig. 
2(e) show histograms of the FT of the calculation and 
the experimental data respectively, the peak positions 
corresponding to the dominant beat frequencies between 




theta (radians) 

FIG. 3: (a) Simulated 'quantum carpet' |ljl] for an ini- 
tially randomly aligned ensemble of room temperature D2 
molecules. For a particular At and 9, the false colour-scale 
indicates the wavepacket probability with a thermal average. 
In the range 230 < At < 290 fs, the ensemble is predicted to 
be maximally aligned with the pump polarization direction, 
(b) The integrated rapid SI yield (normalized to the theory) 
over the same range of At indicates the theory describes our 
results with impressive accuracy: all the major features are 
reproduced. 



rotational states. For a diatomic molecule, AJ — 0, ±2, 
thus the beat frequencies assigned to level J in Fig 2 are 
governed by the spacing between rotational levels J and 
J+2. The peak heights are proportional to the product 
of the respective quantum amplitudes. The nuclear spin 
statistics of D2 are also evidenced by the domination of 
beats between even J states over odd J states. 

The beats for J = - 2, 1 - 3 and 2 - 4 occur at 5.5 THz, 
9.1 THz and 12.8 THz respectively. Rather intriguingly, 
the recent work of Ergler et al 2l| , investigating vibra- 
tional wavepackets in the molecular ion, similar rota- 
tional components have been observed, and attributed to 
the creation of a rotational wavepacket in the molecular 
ion D^. However, the polarizability anisotropy and bond 



length change dramatically under ionization from D2 to 
T>2 , with corresponing beat frequencies being a factor 
of two lower than in D2. The time-domain and equiv- 
alent FT frequency domain data described here are in 
excellent agreement with the expectation of impulsively 
excited rotational wavepackets formed in the neutral D2 
molecule. 

The data presented in Fig. 2 describes only part of 
the rotational dynamics of D2 as those molecules outside 
the angular acceptance of our detector simply are not de- 
tected. To fully understand the emergence of rotational 
order, we turn now to predicting and measuring the rapid 
variation of the angular distribution of molecules around 
the first half-revival. Fig. 3(a) is a false-colour repre- 
sentation of the expected probability distribution as a 
function of evolution time At and 9. When 250 < At 
< 280 fs, maximal alignment is predicted: the ensem- 
ble is expected to be preferentially aligned in a narrow 
double-lobed distribution directed along 9 = 0, tt in this 
'quantum carpet'. By rotating 9 over the range -tt to 
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7r, the axis of revival of the rotational wavepacket is ro- 
tated with respect to our detector axis. The experimen- 
tal data is shown in Fig. 3(b): the relative D + yield as 
a function of 8 and At directly reflects the thermal av- 
erage of the rotational wavepacket probability over the 
thermal average distribution of J. While some statisti- 
cal scatter is present, the features of the quantum carpet 
are well resolved. This demonstrates how the molecular 
axes align through revival to create order. Up to ~ 190 
fs, the ensemble is to a large extent isotropic. However 
in the vicinity of the half revival, the rephasing of the 
wavepacket traps or confines the molecules in a well de- 
fined spatial distribution for a short period of time. It 
should be noted that the large energy spacing of the J- 
states in D2 do not necessitate supersonic cooling, thus 
this demonstration of an ordered molecular ensemble il- 
lustrates a readily accessible quantum system. 

The mapping of this ultrafast rotational wavepacket 
has a number of implications. On the most fundamental 
level, the observation of the quantum carpet in D2 (con- 
firmed by our theoretical predictions) illustrates that a 
room temperature ensemble of this theoretically impor- 
tant molecule can be organized into an ordered state. 
This can allow time- and energy-resolved measurements 
of the electronic and nuclear motion with spectroscopic 
accuracy. Considerable efforts are underway to perfect 
ultrashort (< 5 femtoseconds Fpulse generation through 
molecular phase modulation [10]. By impulsively align- 
ing a high pressure Raman-active gas (short rotational 
period makes D2 a natural choice) in a hollow waveg- 
uide with an ultrafast pump pulse, a locally-aligned ro- 
tational revival propagates along the waveguide at the 
speed of light in the medium, and the molecular motion 
causes a strong phase modulation due to the refractive 
index change. A second pulse (directly equivalent to our 
probe) injected into the waveguide to coincide with the 
rotational revival is then spectrally broadened and com- 
pressed through the phase modulation. Importantly, it 
has been demonstrated that the time dependent phase 
introduced can be accurately controlled by changing the 
delay between the pump and probe pulses. 
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